Effects of interdot hopping and Coulomb blockade on the thermoelectric properties of 

serially coupled quantum dots 
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We have theoretically studied the thermoelectric properties of serially coupled quantum dots 
(SCQDs) embedded in an insulator connected to metallic electrodes. In the framework of Keldysh 
Green's function technique, the Landauer formula of transmission factor is obtained by using the 
equation of motion method. Based on such analytical expressions of charge and heat currents, we 
calculate the electrical conductance, Seebeck coefficient, electron thermal conductance and figure of 
merit (ZT) of SCQDs in the linear response regime. The effects of interdot hopping and electron 
Coulomb interactions on ZT are analyzed. We demonstrate that ZT is not a monotonic increasing 
function of interdot electron hopping strength (t c ). We also show that in the absence of phonon 
thermal conductance, SCQD can reach the Carnot efficiency as t c approaches zero. 
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Introduction 

Recently, many considerable studies have been devoted 
to seeking efficient thermoelectric materials with the fig- 
ure of merit (ZT) larger than 3 because there exist po- 
tential applications of solid state thermal devices such as 
coolers and power generators. [1-6] Some theoretical ef- 
forts have pointed out that a single quantum dot (QD) 
junction system can have a very impressive ZT in the 
absence of phonon conductance. [7-9] However, in prac- 
tice it is difficult to maintain a large temperature gradi- 
ent needed to produce sufficient temperature difference 
across the nanoscale junction. To reduce the tempera- 
ture gradient across the QD junction, it is essential to 
consider many serially coupled QDs.[l,5] The transport 
property of a junction involving N serially coupled QDs 
with strong electron Coulomb interactions is one of the 
most challenging topics of condensed matter physics. To 
gain some insight, we investigate in the present paper 
the thermoelectric effect of serially coupled quantum dots 
(SCQD) as shown in the inset of Figure 1(a). 

It has been shown that the transport properties of the 
SCQD system exhibit several interesting features, includ- 
ing current rectification (due to the Pauli spin blockade), 
negative differential conductance, nonthermal broaden- 
ing of tunneling current, and coherent tunneling in the 
Coulomb blockade regime. [10] Although many theoreti- 
cal investigations of the above phenomena have been re- 
ported, most of them did not investigate the thermoelec- 
tric properties of SCQDs. [11-13] This study investigates 
the ZT of SCQD embedded in a semiconductor nanowire 
with small phonon thermal conductance. [4] It is expected 
that SCQD system has a potential to enhance the ZT 
of nanowires. Here wc consider nanoscale semiconduc- 
tor QDs, in which the energy level separations are much 
larger than their on-site Coulomb interactions and ther- 
mal energies. Thus, only one energy level for each quan- 
tum dot needs to be considered. A two-level Anderson 
model[13] is employed to simulate the SCQD junction 
system. 

Theoretical model 



Using the Keldysh-Green's function technique, [13] the 
charge and heat currents of SCQD connected to metallic 
electrodes are given by 



J 



2c 



deT(e)[f L (e) -f R (e)], 



(1) 



Q = U deT(e)(e-E F -eAV)[f L (e)-f R (e)] (2) 

where T(e) = (7i2(e) + 72i(e))/2 is the transmission fac- 
tor. /i=i(fl= 2 )(e) = l/[e { - t -^^/ kBT ^ R )+l] denotes the 
Fermi distribution function for the left (right) electrode. 
The left (right) chemical potential is given by hl(^r)- 
Tl(r) denotes the equilibrium temperature of the left 
(right) electrode, e and h denote the electron charge 
and Planck's constant, respectively. Te.j(e) denotes the 
transmission function, which can be calculated by eval- 
uating the on-site retarded Green's function (GF) and 
lesser Green's GF.[13] The indices I and j denote the 
QD and the jth QD, respectively. Based on the equation 
of motion method, we can obtain analytical expressions 
of all GFs in the Coulomb blockade regime. Details are 
provided in Ref. 13. The transmission function in the 
weak interdot limit ( t c /U£ <C 1, where t c and Ut de- 
note the electron interdot hopping strength and on-site 
Coulomb interaction, respectively) can be recast into the 
following form 

7Jj(e) = -2 £ ^^Imq im!ff (e), (3) 

where Im means taking the imaginary part of the function 
that follows, and 



G r L m .a( e ) =Vml(H ~ S m ). 



(4) 



t=L(l) 



.i?.(2)( e ) denotes the tunnel rate from the left elec- 
trode to dot A (Ei) and the right electrode to dot B (£2), 
which is assumed to be energy- and bias-independent for 
simplicity. fj,£ = e — En + iTi/2. We can assign the fol- 
lowing physical meaning to Eq. (3). The sum in Eq. 



2 



(3) is over 8 possible configurations labeled by to. We 
consider an electron (of spin a) entering level t, which 
can be either occupied (with probability Nt t „) or empty 
(with probability 1 — Ng_ s ). For each case, the electron 
can hop to level j, which can be empty (with probability 
aj = 1 — Nj >c , — Nj^g + Cj), singly occupied in a spin a 
state (with probability bj & = Nj^ — Cj) or spin a state 



given by 



(with probability bj 



Cj), or a double-occupied 



state (with probability Cj). Thus, the probability fac- 
tors associated with the 8 configurations appearing in 
Eq. (4) become pi = (1 - N e ^)a j7 p 2 = (1 - N l ^ s )b j . s , 

P3 = (1 - N it a)bj t o; Pi = (1 - Nl^)Cj, p 5 = N^dj, 

P6 = Nt t abj t a, Pi = Nt,crbj !<r , and p s = Nt^Cj. S m in the 
denominator of Eq. (4) denotes the self-energy correction 
due to Coulomb interactions and coupling with level j 
(which couples with the other electrode) in configuration 
to. We have Si = t 2 c /fij, S 2 = U Lj + t 2 c /(fij - Uj), E 3 = 
Ue.j+tl/^ij-Uj.i), S 4 = 2U itj +4/(fx j -U j -U j>i ), S 5 = 
U t + ft/fa - U Jjt ), E 6 = U e + U Lj + ft/fa - Uj - U M ), 
Y, 7 = U t + U t j + t 2 J(fj,j - 2Uj, t ), and E 8 = U t + 2U £>j + 
t 2 /(fij — Uj—2Uj y e). E(, Ue, and Uej denote, respectively, 
the energy levels of dots, intradot Coulomb interactions, 
and interdot Coulomb interactions. Here V™ = — 2ImSj 
denotes the effective tunneling rate from level / to the 
other electrode through level j in configuration to. For 
example, T] = -2lmb 2 J ^ = t 2 c Tj/[{e - E.jf + (Tj/2) 2 }. 
It is noted that TJ 1 has a numerator Tj for all configu- 
rations. Furthermore, G r e<J (e) = J2 m =l G £, m .cr( e ) is just 
the on-site single-particle retarded GF for level I as given 
in Eq. (A16) of Ref. 13, and G\ ma {e) corresponds to its 
partial GF in configuration to. The transmission func- 
tion written this way has the same form as Landauer's 
formula for a single QD with multiple energy levels in- 
cluding intralevel and intcrlevel electron Coulomb inter- 
actions. [14,15] 

The probability factors of Eq. (3) are determined by 
the thermally averaged one-particle occupation number 
and two-particle correlation functions, which can be ob- 
tained by solving the on-site lesser Green's functions, [13] 



and 

ci = 
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de * T e f e (e)+Tffj(e) 
7T ^ Te + TT 



JmGl m Je). (6) 



Note that I ^ j in Eqs. (3), (5), and (6). In the linear 
response regime, Eqs. (1) and (2) can be rewritten as 



J - L n-jT +^12-^2 

AV AT 
Q = £-21^7- + £22^2", 



(7) 

(8) 



where AV = ^l— A* R an( i AT = Tj, — Tr are the voltage 
and temperature differences across the junction. Ther- 
moelectric response functions in Eqs. (7) and (8) are 



r 2e 2 T f df(e) 

Cii = —J deT{e) W )t ' 



2eT 2 



deT(e)( 



9_m 

dT 



(9) 



(10) 



^i = ^fj deT(e)(e - E F )(^) T , (11) 



and 



£22 = — J deT{e)(e - E F )(-^-)i 



(12) 



Here T(e) and /(e) = l/[ e ^-E F )/k B T + ^ are evaluated 
in the equilibrium condition. It can be shown that the 
Onsager relation C12 = £21 is preserved. These thermo- 
electric response functions can also be found in ref. [7], 
where authors investigated the thermoelectric properties 
of a single QD. 

If the system is in an open circuit, the electrochemical 
potential will form in response to a temperature gradi- 
ent; this electrochemical potential is known as the See- 
beck voltage (Seebeck effect). The Seebeck coefficient 
(amount of voltage generated per unit temperature gra- 
dient) is defined as S = AV/AT = -£ 12 /(T£ n ). To 
judge whether the system is able to generate power or 
refrigerate efficiently, we need to consider the figure of 
merit, 1 which is given by 



ZT = 



S 2 G e T _ (ZT) C 



~h ftph 1 H~ Kphf K>e 



(13) 



Here G e = Cu/T is the electrical conductance and 
K e = {{£22 /T 2 ) — £nS 2 ) is the electron thermal con- 
ductance. (ZT)o represents the ZT value in the absence 
of phonon thermal conductance, n p h- For simplicity, we 



Kph = KphfiFs- [16-18] K p h,0 



3h 



is the uni- 
versal phonon thermal conductance arising from acoustic 
phonon confinement in a nanowire, [16-18] which was con- 
firmed in the phonon wave guide. [19] The expression of 
Uph = KphflF s with F s = 0.1 can explain well the phonon 
thermal conductance of silicon nanowire with surface 
states calculated by the first-principles method. [16] The 
dimensionless scattering factor F s arises from phonon 
scattering with surface impurities or surface defects of 
quantum dots. [1,16] Here, we adopt ^.=0.02, which is 
smaller than F s =0.1, because QDs can enhance the 
phonon scattering rates and reduce phonon thermal con- 
duction as pointed out in Ref. 1. 
Results and discussion 

Here, we consider the case of identical QDs in the op- 
timization of ZT, although it is understood that the size 
fluctuation of QDs can suppress .ZT. [13] In Fig. 1(a) and 
1(b), we plot (ZT)o and ZT as a function of temperature 
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for various electron hopping strengths. We adopt the fol- 
lowing physical parameters: Eg = Ep + 30r , Ug = 30r , 
Uej = 10r , and T L = T R = V = IT . All energy scales 
are in the units of the characteristic energy, Tq. In Fig. 
1(a), we see that (ZT)o increases with decreasing t c and 
diverges as t c — > 0. This behavior can be proved rigor- 
ously as we shall illustrate below. It implies that SCQD 
can reach the Carnot efficiency in the limit of extremely 
weak interdot coupling, if one can fully suppress K p h, for 
example by inserting a nanoscale vacuum layer to block 
the phonon heat current. Although it would be a chal- 
lenging task to implement a vacuum layer between one of 
the electrodes and SCQD, it may be possible to test this 
idea out via a scanning tunneling microscopic (STM) ex- 
periment by using a set-up as shown in the inset of Fig. 
1(a). In Fig. 1(b), we see that ZT is enhanced with in- 
creasing t c until t c reaches 3Fo, and it becomes reduced 
for higher t c . 

The diverging behavior of (ZT)q with respect to t c 
is further illustrated by the results of Fig. 2(d). The 
maximum ZT is suppressed in the presence of n p h, which 
is much larger than n e for small t c . The behaviors of ZT 
shown in Fig. 1(b) are mostly determined by the power 
factor (S 2 G e ). Once t c is larger than 3Fo, the reduction 
of S 2 is faster than the increase of G e . This explains why 
the maximum ZT at t c = 4Tq is smaller than that at t c = 
STq. The location of ZT max is nearly independent of t c , 
and it occurs near k B T = 8.8r . For comparison, we also 
show the results (curves with triangle marks) for the case 
without electron Coulomb interactions in Fig. 1(b). It 
is seen that the maximum ZT is enhanced when we turn 
off the electron Coulomb interactions. Such a behavior 
is similar to that of a single QD with multiple energy 
levels. [7,8] The effect of electron Coulomb interactions is 
significant only for temperature between 6Fo and 50Fo. 
Namely, the electron Coulomb interactions are negligible 
when U/(k B T) > 1 or U/(k B T) < 1. 

To further understand the behavior of ZT with respect 
to t c , we plot the electrical conductance (G e ), Seebeck 
coefficient (S 1 ), electrical conductance K e , and (ZT)q as 
functions of t c in Fig. 2 for various detuning energies, 
A = Ei — Ep. When Eg is close to the Fermi energy, 
G e and K e are enhanced, whereas S and (ZT)q are sup- 
pressed. The behavior of (ZT)q at A = 30Fo in the ab- 
sence of Coulomb interactions is also shown by the curve 
with triangles, which has similar trend as the solid line. 
Thus, it is instructive to analyze (ZT)q in the absence 
of Coulomb interactions. Keeping the leading order of 

t 2 , we have £ u = jj^^fc cos h=(A/2fc B T)> C ^ = £ 2i = 

2e *c A j r _ _J2__ti A 2 

hk B T /2 cosh 2 (A/2k B T) ' mlLl / " 22 — hk B To/2 cosh 2 (A/2k B T) • 

Therefore, G e oc t 2 , S = —k B A/eT is independent on 
t c , and k e — (£ 2 2 — C\ 2 / Cn)/T 2 vanishes up to t 2 . 
Thus, the leading order of K e is t\. This indicates that 
(ZT)q oc 1/t 2 in the limit of weak interdot hopping. 

Figure 3 shows ZT as a function of A = Eg — Ep for 
various electron hoping strengths at fc^T = lOFo. Other 
physical parameters are kept the same as those for Fig. 



1. When t c = O.lTo, the maximum ZT (ZT max ) occurs 
at near A = 27Fo. The peak position only shifts slightly 
to higher A with increasing t c . We have ZT max = 2.79 
and 3.18 for t c = ITq and 3Fo, respectively. However, at 
t c = ATq we have ZT max = 3.07, which is smaller than 
ZT max for t c = 3Tq . Thus, it also illustrates that ZT is 
not a monotonically increasing function of t c . We further 
calculated ZT as a function of t c for A = 10, 20, 30Fo and 
k B T = IOTq in the presence of k p h and found that again 
ZT is not monotonically increasing function of t c (not 
shown here) . We conclude that as long as K p h dominates 
over K e , the t c dependence of ZT is mainly determined 
by the power factor S 2 G e , where the behaviors of G e 
and S are similar to the results shown in Fig. 2(a) and 
2(b). When i c /To < 1, G e increases much faster that the 
reduction of S 2 for increasing t c , and the power factor 
slowly reaches the maximum when t c approaches 3Fo. 
When t c > 3Tq, the power factor decreases due to the 
fast reduction of S 2 which prevails over the increase of 
G e . The curve with triangle marks is for t c = 3Tq in 
the absence of Coulomb interaction. We see that ZT max 
is larger when Ug = Ui.j = 0. Based on the results of 
Fig. 3, we conclude that it is important to control the 
detuning energy A for the optimization of ZT . 

In Figures 1,2 and 3 we have considered the case with 
Ep below QD energy levels. It would be interesting to 
investigate the case with Ep above the energy levels of 
QDs. Fig. 4 shows G e , S, K e , and ZT of an SCQD 
with t c = 3Fo as functions of applied gate voltage for 
various temperatures. Once t c > (Fl + Fr) = 2Fo, the 
eight peaks for G e can be resolved at k B T = lTo. These 
eight peaks correspond to the following resonant chan- 
nels: Et-t c , E e + t c , E t + Uej -t c , E e + Utj + t c , 
Ep+Ui^+Ue-tc, Ei+Ue^+Up+tc, E e +2U ej +U e -t c , and 
Ei+2Uc.j + U(+t c , which are tuned by the gate voltage to 
be aligned with Ep. These eight channels result from the 
four configurations of Pi, Ps, P6, and p$ in Eq. (4). Such 
a result implies that SCQD with identical QDs acts as a 
QD with effective two levels of Eg — t c and Eg+t c and sat- 
isfying Hund's rule. These eight peaks are smeared out 
with increasing temperature. The sign changes of S with 
respect to the gate voltage result from the bipolar effect, 
i.e., the competition between electrons and holes, where 
holes are defined as the unoccupied states below _Ep\[13] 
The electronic thermal conductance (n e ) also exhibits 
eight peaks, and we noticed that the local maxima of the 
K e curve nearly coincide with the local minima of the G e 
curve. We see that ZT values are still larger than 3 even 
when Eg is deeply below Ep (say, at eV 9 = 70Fo). This is 
attributed to the electron Coulomb interaction. To illus- 
trate that, we also show the results with Ui = Ugj = at 
k B T = STo (see the curve with triangle marks). The os- 
cillation of ZT in the case of Up = Ugj = is attributed 
to the sign change of S at V g = 10T . Note that S goes 
to zero at V g = 10T , which results from the electron- 
hole symmetry (with Eg + t c and Eg — t c straddling Ep 
symmetrically). We see that ZT vanishes for eV g > 40r 
in the absence of electron Coulomb interactions. Unlike 
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the case of Ep < Eg, where the finite U causes reduction 
of ZT, here the electron Coulomb interaction leads to 
enhancement of ZT when Ep > E( . 
Conclusions 

In summary, the thermoelectric properties including 
G e , S, K e and ZT of the SCQD junction system are in- 
vestigated theoretically We demonstrate that the Carnot 
efficiency can be reached when t c approaches zero in 
the absence of phonon thermal conductance. When the 
phonon contribution dominates the thermal conductance 
of SCQD junction, the optimization of ZT can be ob- 



tained by the thermal power defined as S 2 G e . We also 
found that the presence of electron Coulomb interactions 
can lead to either reduction or enhancement of ZT de- 
pending on whether the Fermi level is below or above the 
QD level. 
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Figure Captions 

Fig. 1. (a) (ZT) and (b) ZT as functions of tem- 
perature for various interdot hopping strengths (t c = 
0.1,0.5,1,3, and 4r ). E e = Ep + 30r , U e = 30r , 

u td = ior , and r L = r = r . 

Fig. 2. Electrical conductance (G e ), Sccbeck coeffi- 
cient (S), electrical thermal conductance (n e ), and (ZT)o 
as functions of t c at fcgT = 5r for A = 10r (dotted 
curves), 20Fo (dashed curves), and 30Fo (solid curves). 
Other parameters are the same as those of Fig. 1. 

Fig. 3. ZT as a function of A for different electron 
hopping strength at ksT = lOFo. Other parameters are 
the same as those of Fig. 1. 

Fig. 4. G e , S, K e , and ZT as a function of applied 
gate voltage for fc^T = ir (solid), 2r (dashed), and 
3r (dotted). E t = E F + 10r o and t c = 3T a . Other pa- 
rameters are the same as those of Fig. 1. The curves with 
triangle marks are for the case without electron Coulomb 
interactions for fc^T = 3Fq . 
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